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ABSTRACT 

This paper was to verify the STAGS (general shell, geometric and material nonlinear) 
code and the critical crack-tip-opening angle (CTOA) fracture criterion for predicting 
stable tearing in cracked panels that fail with severe out-of-plane buckling. Materials 
considered ranged from brittle to ductile behavior. Test data used in this study are 
reported elsewhere. The STAGS code was used to model stable tearing using a critical 
CTOA value that was determined from a cracked panel that was “restrained” from 
buckling. The analysis methodology was then used to predict the influence of buckling on 
stable tearing and failure loads. Parameters like crack-length to-specimen-width ratio, 
crack configuration, thickness, and material tensile properties had a significant influence on 
the buckling behavior of cracked thin-sheet materials. Experimental and predicted results 
showed a varied buckling response for different crack-length-to-sheet-thickness ratios 
because different buckling modes were activated. Effects of material tensile properties and 
fracture toughness on buckling response were presented. The STAGS code and the CTOA 
fracture criterion were able to predict the influence of buckling on stable tearing behavior 
and failure loads on a variety of materials and crack configurations. 

NOMENCLATURE 

B Specimen thickness, mm 

q Initial half-length of crack, mm 

d Minimum element size along crack line, mm 

E Young’s modulus, GPa 

* 
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Fb Buckling response ratio (Pb/Pnob) 

Pb Peak load when sheet allowed to buckle, N 

Pnob Peak load when sheet restrained from buckling, N 

S Applied stress, MPa 

Sk Peak stress during crack linkup (k = 1 to 3), MPa 
Sf Failure stress, MPa 

W Out-of-plane displacement in z-direction, mm 

w Specimen half-width, mm 

x,y,z Cartesian coordinates 

Ac Crack extension, mm 

V Poisson’s ratio 

o X x Normal stress in x-direction 

y c Critical crack-tip opening angle, degrees 

Critical opening angle for sawcut at crack initiation, degrees 

INTRODUCTION 

Stably tearing cracks in fuselage structures are subjected to complex loading 
conditions (pressure, biaxial, and shear) and cracks in fuselages are subjected to mixed- 
mode and out-of-plane deformations. In the framework of the NASA Airframe Structural 
Integrity Program, the FRANC3D/STAGS codes are being developed to predict the 
behavior of these cracks under complex fuselage loading with the critical crack-tip- 
opening angle (CTOA) fracture criterion [1]. In developing the residual strength methods, 
the ability of STAGS to predict the influence of bulging or buckling on stable tearing 
behavior, in the presence of multiple-site and discrete source damage, needed to be 
verified on laboratory tests and on structural test articles. Development of these residual 
strength analysis methodologies for complex loading are essential to ensure the continued 
safe operation of the aging commercial transport fleet and aircraft for the next generation. 
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Several fracture criteria have been proposed to predict the residual strength of 
cracked specimens or structures using the finite-element method. The crack-tip stress or 
strain, crack-tip-opening angle (CTO A), crack-tip-opening displacement (CTOD), crack- 
tip force, strain-energy-release rate, and the J-integral have been employed to predict the 
stable crack growth behavior under monotonic loading. Of these, the CTOA or CTOD 
criterion had been shown to be well suited for modeling the stable crack growth process 
and instability during fracture [2-4]. Experimental measurements of the critical CTOA 
during the fracture process for thin-sheet aluminum alloys [5] have shown that CTOA 
remains nearly constant during stable crack growth after a small amount of crack 
extension (about one thickness). 

Because cracks in fuselage structures are subjected to out-of-plane deformations, the 
influence of bulging or buckling on stable tearing must be accounted for in the analysis. 

By addressing a simpler crack problem, such as the fracture behavior of middle-crack 
tension M(T) specimens under buckling conditions, an assessment of the capability of the 
STAGS code and the critical CTOA criterion can be made. 

In 1962, Kuhn and Figge [6] developed an expression to determine the buckling 

response ratio Fb (ratio of buckling load to failure load of restrained M(T) specimens with 

anti-buckling guides) as a function of the initial crack-length-to-thickness (q/B) ratio. 

Later, Forman [7] carried out a series of experiments to determine the effects of buckling 
and specimen dimensions on residual strength of thin- sheet aluminum alloys and steels. At 
least two M(T) specimens were tested with guide plates lightly clamped against the panels 
to prevent buckling and one specimen was tested without guide plates to allow buckling to 
occur. The test program showed that there was a considerable reduction in residual 
strength of cracked sheets which were allowed to buckle. From these results the linear 
relationship given by Kuhn and Figge, to define the buckling response ratio, could not be 
generalized for most materials. This lack of uniqueness in representing the buckling 
response ratio with thickness for different materials, and the need to verify the STAGS 
code and the CTOA fracture criterion, prompted this study. 
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The objective of this study was to analyze cracked panels that fail with severe out-of- 
plane buckling using STAGS (STructural Analysis of General Shells) code [8,9] and the 
CTOA fracture criterion. Materials with fracture behavior ranging from brittle to ductile 
were considered. The STAGS code was used to model stable tearing using a critical 
CTOA value that was determined from a cracked panel that was “restrained” from 
buckling. The STAGS code and the critical CTOA value were then used to predict the 
influence of buckling on failure loads. The influence of crack-length-to- specimen- width 
ratio, thickness, crack configuration and material tensile properties on the buckling 
behavior on several thin-sheet aluminum alloys and a steel were studied. These materials 
covered a wide range in fracture toughness behavior and had critical CTOA values ranging 
from a very low value of 1 .9 degrees to a value of 4.6 degrees. Experimental test data used 
herein are reported elsewhere. A hypothetical material with a 10 degree critical CTOA 
value was also considered to study the buckling behavior of specimens made of an 
extremely high toughness material. Comparisons are made between the measured and 
predicted applied stress against crack extension, in-plane crack-opening displacement 
(COD) and crack out-of-plane displacement (COPD) for 2024-T3 specimens tested with 
and without buckling guides. The influence of various buckling modes and the initial 

crack-length-to-thickness ratio (q/B) on failure loads were studied. The effects of material 

tensile properties and fracture toughness on buckling response are presented and discussed. 
The influence of multiple-site-damage (MSD) on the buckling response were also studied 
for some previously conducted tests on 2024-T3 panels. (MSD is the presence of a large 
number of small fatigue cracks near a large primary lead crack, such as fatigue cracks along 
a rivet row in a lap-splice joint in a fuselage.) The ability of the STAGS shell code and the 
CTOA fracture criterion to predict the influence of buckling on stable tearing behavior and 
residual strength of thin-sheet materials were investigated. 

MATERIALS AND SPECIMEN CONFIGURATIONS 

Analyses were conducted on middle-crack tension M(T) specimens and on multiple- 
site damage (MSD) specimens subjected to remote tension. The materials studied were 
2024-T3, 2219-T87 and 7075-T6 aluminum alloys and an AM355CRT steel. These 
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materials cover a wide range in fracture behavior from brittle to ductile. In general, the 
initial crack-length-to-half- width (ci/w) ratio of 1/3 was considered but, in addition, the 

2219-T87 material was analyzed with Cf/w ratios of 0.2, 0.4, 0.6 and 0.8 to study the 

effects of crack length on residual strength. As shown in Table 1, the materials stress- 
strain behavior were approximated by piece- wise-linear curves. 

FINITE ELEMENT ANALYSES 

In the M(T) specimen, the normal stress, g xx parallel to the crack surface is 

compressive, as shown in Figure 1. If the specimen is not restrained, then the crack flanks 
will deform out-of-plane (buckle) and this behavior will greatly reduce the load carrying 
capacity. The STAGS code and the CTOA fracture criterion will be used to study stable 
tearing and buckling interactions in this specimen type. 


s 



Figure 1. Middle-crack tension, M(T), specimen. 
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STAGS is a finite element program for the analysis of general shell-type structures [8, 
9]. The program has several types of analysis capabilities (static, dynamic, buckling, crack 
extension, material nonlinear and geometric nonlinear behavior). STAGS has crack 
extension capability based on the critical crack-tip-opening angle or displacement (CTOA 

or CTOD) criterion and the traditional KR-curve. Its finite element library includes, two- 

node beam elements, triangular shell, quadrilateral shell and transitional elements to 
connect triangular and quadrilateral elements (see refs. 8 and 9 for more details). In the 
current study, quadrilateral shell elements with 6 degrees-of-freedom per node (three 
displacements and three rotations) were used. The quadrilateral shell element was under 
“plane-stress” conditions. Elastic-plastic material behavior was approximated by multi- 
linear stress-strain curves (see Table 1). The White-Besseling plasticity theory with a form 
of kinematic hardening was used to account for yielding and reverse yielding during 
unloading [9]. 

The specimen configuration and typical mesh pattern are shown in Figures 1 and 2, 
respectively. Because the configuration and loading were symmetric, only a quarter of the 
sheet was modeled, as shown in Figure 2. The minimum crack-tip element size (d) of 1 
mm, along the line of crack extension, was chosen to be the same for all meshes generated 
for different case studies. From previous parametric and convergence studies, it was 
found that a minimum crack-tip element size of 1 mm (linear-strain element) was sufficient 
to model stable tearing under elastic-plastic conditions [10]. Crack growth was governed 

by monitoring the critical CTOA (vj/ c ) at one element distance (1 mm) behind the crack tip. 

In general, the critical value (x\f c ) was determined by matching the average failure load 

measured on several tests of M(T) specimens (restrained from buckling), as will be 
discussed later. 

A critical CTOA (\|/ c ) was used to evaluate the onset of crack growth. This criterion 
is equivalent to a critical CTOD (8 C ) value at a specified distance behind the crack tip. At 

each load increment, the CTOA was calculated and compared to a critical value \[/ c . 

When the CTOA exceeded the critical value, the crack-tip node was released and the 
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crack was advanced to the next node. This process was continued until crack growth 
became unstable under load control or until the desired crack length had been reached 
under displacement control. 


Y 

Z X 



Figure 2. A typical finite-element model of one-quarter of an M(T) specimen. 

The aerospace industry frequently uses sawcuts to approximate fatigue cracks 
because they are less costly to manufacture and they can be placed more accurately in 
complex structural test articles than an actual fatigue crack. To approximate the use of 
sawcuts, the sawcut was assumed to undergo a deformation at the sawcut tip before a 

crack would initiate. In order to simulate the sawcut effects, a critical sawcut angle, \|/j , 

was selected to model the deformations that take place at the sawcut tip before a crack 
would initiate. Again this value was selected to match the load required to initiate a crack 
at the sawcut. In the analysis, once the critical sawcut angle was reached, the node at the 
tip of the sawcut was released and the crack advanced to the adjacent node (1 mm 
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spacing). Thus, the crack extension at initiation was 1 mm. After a crack had initiated, 
the critical tearing CTO A, y c , was used to grow the crack. The tearing angle, \|/ c , and 

sawcut angle, \|/j , for the different materials considered in the present analyses are listed 
in Table 2. 

In order to simulate buckling in the STAGS numerical analyses, a concentrated load 
of a very small magnitude was gradually applied in small steps in the transverse or out-of- 
plane (z) direction. The concentrated load was applied at the center of the crack (x = 0). 
A trial-and-error procedure was used to find the small load that was sufficient to induce 
buckling. The out-of-plane load was gradually reduced as the external loading was 
applied. 

RESULTS AND DISCUSSIONS 

In the following sections, the results of the numerical analyses are compared with the 
experimental observations and discussed in detail. To begin with, stable crack growth test 
results, measured crack opening displacements (COD) and measured crack out-of-plane 
displacements (COPD) are compared with numerical results from fracture simulations. 
Later, the sawcut simulation capability in STAGS is discussed and analyses are compared 
with some test results. Finally, the effects of buckling on stable tearing and the influence 
of multiple-site damage (MSD) cracking on crack growth behavior with buckling are 
presented. 

Fracture Simulations with Fatigue Cracks 

The numerical analyses of M(T) specimens with and without buckling and the 
corresponding test data are plotted as applied stress against crack extension in Figure 3. 
Square and circular symbols represent the tests conducted with guides to prevent buckling 
and the triangles show data where guides were not used and the sheet was allowed to 

buckle [11]. An angle of \jr c = 4.6 degrees was used to fit the average of the experimental 

failure loads corresponding to the specimens with anti-buckling guides (restrained). The 
dashed curve corresponds to the fracture simulation of the specimens where out-of-plane 
displacements were restrained. Later, the analysis was repeated with the same angle but 
buckling was allowed. The solid curve represents the numerical results when out-of-plane 
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deformations were permitted (unrestrained). Both the tests and numerical results indicate 
a significant influence of buckling on stable tearing. The residual strength was drastically 
reduced when the sheet was allowed to buckle. The experimental and numerical results 
agreed well near the maximum applied stress but both analyses predicted more crack 
extension in the early stages of crack growth than in the tests. This discrepancy has been 
traced to severe crack tunneling [4], since crack extension was measured at the free 
surface, and to the plane-stress behavior assumed in the analysis [5, 12]. However, these 
results demonstrate that STAGS and the CTOA criterion can model the influence of 
buckling on failure loads for cracks in thin-sheet aluminum alloys. 



Figure 3. Measured and predicted stable crack growth in M(T) specimens using CTOA. 

Comparison of Measured and Predicted COD — A comparison of applied stress against 
crack-opening displacement (COD) from the test (symbols) and the STAGS analysis is 
shown in Figure 4. The crack-opening displacements in the y-direction were measured at 
the centerline of the crack. Only test results with anti-buckling guides were available for 
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comparison [11]. Here, there was good agreement between the tests and numerical results 
(dashed curve) for the restrained specimens over the complete test range. The solid curve 
shows the STAGS numerical analyses when out-of-plane displacements were allowed. 
When out-of -plane deformations were allowed (buckling), the calculated COD values 
were much larger for a given applied stress than the restrained specimens. 

Measured and Predicted Out-of-Plane Displacements — The crack out-of-plane (or 
transverse) displacements (COPD), measured at the center of the crack (x = 0) in the z- 
direction, are plotted against applied stress in Figure 5. The symbols show the test 
measurements [11] and the solid curve shows the results from the numerical analysis. 

After reaching the maximum stress, the numerical results predicted larger COPD values 
than the test results as the specimen unloaded. This behavior may be attributed to 
constraint effects because the state-of-stress in the STAGS analysis was plane stress. 
Higher constraint conditions, such as plane strain, causes a more rapid reduction in the in- 
plane displacements after maximum load [12]. Thus, it was expected that the COPD 
trends would exhibit a similar behavior after maximum load. Overall the results from the 
analysis showed good agreement with the test data up to the maximum stress. 



Figure 4. Measured and predicted crack-opening displacement (COD) for M(T) specimen. 
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Figure 5. Measured and predicted crack out-of-plane displacement (COPD) 
for M(T) specimen. 

Fracture Simulations with Sawcuts 

Comparison of experimental and numerical applied stress against crack extension for 
an M(T) specimen with a sawcut are shown in Figure 6. The test results (symbols) are 
from a restrained M(T) specimen with a sawcut (1.8 mm wide sawcut blade). As 
previously mentioned, the behavior of sawcuts are quite different from that of a fatigue 
crack. The sawcut had a blunt tip and required much higher applied stresses to initiate a 

tearing crack. In order to simulate sawcut effects, a critical sawcut angle, \| q , was 

selected to match the applied stress required to initiate a crack at the tip of the sawcut, 
similar to that developed in reference 10. Once a crack had initiated, the critical tearing 

angle, \j/ c , was used to grow the crack. 

Several types of fracture analyses were performed with STAGS. First, a tearing crack 
was analyzed with a critical angle of 4.6 degrees (dash curve). Then analyses were 

conducted to find, by trial- and-error, the critical sawcut angle. A sawcut angle of \}/i = 15 
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degrees was found from these analyses. The solid curve shows the STAGS analysis with 
the sawcut angle to initiate a crack and a critical angle of \|/ c = 4.6 degrees to tear the 

crack. The fit to the experimental results was reasonable. As was evident from Figure 6, 
the maximum applied stress was not greatly affected by having either a sawcut or crack. 
But there was a significant difference (about 100 %) in the initial applied stress required to 
initiate a crack at the sawcut and the applied stress required to initiate a tearing fatigue 
crack (see curves at Ac = 0 in the Fig. 6). 



Figure 6. Comparison of measured and calculated stable crack growth 
in M(T) specimen with sawcut or crack using CTOA. 

Effects of Buckling on Stable Tearing 

Earlier experimental studies [6, 7] have shown that parameters like crack-length-to- 
width ratio (q/w), thickness (B), and material tensile properties have a significant 

influence on the buckling behavior of M(T) specimens. In order to understand the 
influence of these parameters, a wide range of materials were considered in the present 
study (2024-T3, 2219-T87, 7075-T6 aluminum alloys and AM355CRT steel). These 
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materials have critical CTOA values from 1 .9 to 4.6 degrees. A hypothetical material with 
a 10 degree critical angle was also considered in the analyses to study the buckling 
behavior of specimens made of a very ductile material. 

Residual Strength of Ductile Materials — In Figure 7, the buckling ratio (Fb = Pt/Pnob) 
from both experiments and numerical analyses are plotted against the initial crack length- 
to-thickness (q/B) ratios. Fb is the ratio of the peak load (Pb) when the thin sheet was 

allowed to buckle to the load (P n ob) when all out-of-plane displacements are restrained. 
The results shown in Figure 7 were tested with specimens that had the same crack-length- 
to- width ratio, q/w =1/3. The specimen widths (2w) ranged from 75 to 305 mm for the 
1.5-mm thick aluminum alloy and steels [7] and ranged from 75 to 610 mm for the 2.3- 
mm thick aluminum alloy [11]. For larger values of the q/B ratio, the buckling ratio 

decreased and tended to level off at about 0.7. Both the aluminum alloys and steels 
exhibited nearly the same behavior. 



Figure 7. Measured and predicted results for 2024-T3 aluminum alloy and 
AM355CRT steel. 
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For comparison, the linear approximation from Kuhn and Figge [6] was shown as the 
dashed line. The experimental results from the tests [7, 1 1] for both materials lie on either 
side of the linear approximation. Thus, the linear relationship can only be considered as a 
first approximation to the buckling response. Actually, the experimental results show that 
the behavior was nonlinear. 

The STAGS numerical results are shown by solid and dash-dot curves for the 
aluminum alloy 2024-T3 and AM355CRT steel, respectively. Here the fracture analyses 

were conducted at various q/B ratios by varying B with q/w as constant (q/w = 1/3). 

Numerical results compared well with the experimental trends. Both the aluminum alloy 
and steel predictions show similar buckling response. Their corresponding critical angles 

(Vc) were 4.6 and 4.4, respectively. With an increase in the q/B ratio from 0 to 50, the 

buckling ratio dropped linearly from unity to about 0.8. For q/B ratios greater than 50, 

the buckling ratio decreased slowly and approached a value of about 0.7, but the buckling 
ratio oscillated. This nonlinear behavior was suspected to be due to the buckhng response 
shifting from one buckling mode to another. In the initial stage, the buckling response was 

in the first buckling mode. With an increase in q/B, the sheet reached a critical position 

and the response shifts into the second buckling mode. This change from the first mode to 
a higher buckling mode provide an increased resistance to buckling and causes the ratio to 

decrease at a slower rate. For a q/B ratio of 200, the sheet exhibited severe out-of-plane 

deformations both parallel and perpendicular to the crack surfaces. 

In an effort to study the influence of the shift from one buckling mode to another on 

residual strength, the ST AGS/CTO A analysis was conducted on panels with q/B ratios 

less than 50, greater than 50 and much larger than 50. The initial crack length was the 
same for all analyses. The magnitude and direction of the initial out-of-plane load to 
initiate buckling was the same in all cases. The out-of-plane W-displacement (at the peak 

stress) is plotted against the x/w value in Figure 8. The results at the smallest q/B value 
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shows the first mode of buckling. For larger q/B values, the second and third mode of 
buckling developed. 

Residual Strength of Brittle Materials — Forman [7] also conducted fracture tests on 
various specimens made of 7075-T6 aluminum alloy (2w ranged from 38 to 305 mm) with 

and without buckling guides. The buckling ratio (Pb/Pnob) from these tests are shown as 

symbols in Figure 9 as function q/B. For the small specimen widths tested, the buckling 

ratio was, surprisingly, greater than unity on several test series. The reason for this 
behavior may be specimen scatter because this material is quite brittle in comparison to the 
2024-T3 alloy. Again, the dashed line is the simple buckling equation proposed by Kuhn 
and Figge [6]. Here the Kuhn-Figge expression was in good agreement with the test data. 



Figure 8. Predicted out-of-plane buckling mode shapes along crack line 
for 2024-T3 aluminum alloy M(T) specimens. 


An analysis with the STAGS code on the restrained 7075-T6 specimens with the 
largest width produced a critical angle of 1.9 degrees. When the buckling analyses were 
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repeated with the 7075-T6 alloy, the buckling response was also linear in the range of 0 < 
Cj/B < 100, as shown in Figure 9, but the slope was quite different than that for the 2024- 

T3 alloy and steel, as shown in Figure 7. This material exhibited very small amounts of 
stable tearing and had very little plastic deformation at the crack tip before fracture. 
However, there was good correlation between the test results and the numerical analyses 
with the CTOA concept. These results demonstrate that the buckling response is a strong 
function of the fracture toughness of the material. 



Figure 9. Measured and predicted behavior for 7075-T6 aluminum alloy. 

Residual Strength of Hypothetical ('Very Ductile ) Material — To further study the 
influence of toughness on the residual strength during buckling, a 2024-T3 aluminum alloy 

with a hypothetical critical angle, \|/ c , of 10 degrees was analyzed. These results are 

shown in Figure 10. As previously mentioned, this angle was selected to study the effect 
on buckling when the size of the plastic zone becomes very large compared to the 
specimen width. For comparison, the residual strength for the same material with a critical 
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angle (\|/ c ) of 4.6 degrees is also shown. The analyses showed that buckling had much less 

influence on the failure loads for the very ductile alloy because the material was failing 
under net- section yielding conditions. The drop in the buckling ratio over the wide range 

of q/B ratios analyzed was only about 10 percent. 



Figure 10. Measured and predicted behavior for two critical CTO A values. 

Effect of CTOA on Residual Strength — From Figure 10, it was found that the residual 
strength of the 2024-T3 aluminum alloy was much less when the critical angle \|/ c was 4.6 

degrees instead of 10 degrees. In order to study the tearing-buckling behavior as a 
function of the critical angle, the buckling response of the 2024-T3 alloy was analyzed 

with various critical angles from 2 to 10 degrees with a q/B ratio of 200. These results 

are shown in Figure 11. At high angles, the material fails by net-section yielding and 
buckling had a small influence failure loads (about 10 %). At the lower angles, the 
material had low fracture toughness and failure was brittle in nature and, again, buckling 
had a small influence failure loads. However, for CTOA values between these two 
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extreme cases (3 to 5 degs.), buckling effects on stable tearing and failure loads were 
significant. These analyses indicate that buckling response is a strong function of fracture 
toughness or critical angle. 



Figure 11. Variation in buckling response ratio with CTO A. 

Effect of Crack Length, Width, and Thickness on Residual Strength — In this section, 
the influence of various crack lengths, widths, and thicknesses on buckling response is 
presented and discussed. Earlier, the buckling response for 2024-T3 material was 

obtained at various values of q/B by varying B and keeping q/w as constant (Fig. 7). In 
the present case, the thickness B was held constant (1.9 mm) and both q and w were 
varied while maintaining q/w =1/3. The results from these analyses are presented in 

Figure 12. For low q/B ratios, the buckling response was nearly the same by varying 
either B or w. But the differences in the buckling response are significant in specimens 
with q/B ratios greater than 50. These results suggest that larger panels tend to buckle 

more and produce lower failure loads than smaller width panels at the same q/B ratio. 
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Figure 12. Buckling response ratio in large panels for various thicknesses, 
crack lengths and widths. 

Effect of Crack-Length-to-Width Ratio on Residual Strength — The effects of varying 
the crack-length-to-width ratio, q/w, on buckling response is shown in Figure 13. Three 

510-mm wide (2w) panels with various crack lengths were analyzed with buckling to 
determine the critical angle of 3.6 degrees (unpublished data from Royce Forman, NASA 
Johnson Space Center). The 2219-T87 aluminum alloy was then analyzed with STAGS 

using the critical angle of 3.6 degrees for q/w ratios of 0.2, 0.4, 0.6 and 0.8 with 2w = 

510 mm. There were no corresponding test data for the influence of q/w on buckling 
response. The numerical results show that the buckling response was a slight function of 
q/w with variations of about 10 percent for the various q/w ratios analyzed. Also, the 
analyses indicate that the activation of the higher buckling modes, as shown by the 
oscillatory nature of the buckling response, was weakly dependent upon the q/w ratios. 
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Figure 13. Buckling response ratio as a function of initial crack-length-to- width 
ratio for 2219-T87 material. 

Influence of MSD Cracking on Stable Tearing and Buckling 

Small multiple-site-damage (MSD) cracking at adjacent structural details, such as in a 
riveted fuselage structure, have a strong influence on residual strength, see Swift [14]. In 
the previous sections, the effects of buckling on stable tearing in the presence of a single 
large crack was presented and discussed for various materials. Previously, Newman et al. 
[13,15] studied the influence of MSD on stable crack growth behavior in 2024-T3 panels. 
In order to study the influence of buckling on MSD crack behavior, a STAGS analysis was 
carried out on 2024-T3 MSD panels which had a single large lead 
crack with two-symmetric MSD cracks on either side of the lead crack (crack number 3), 
as shown in Figure 14. 

Broek et al. [16] conducted fracture tests on 508-mm wide M(T) and MSD 
specimens made of 2024-T3 Alclad material. Sawcuts were used instead of fatigue 
cracks. The M(T) specimens with a single sawcut were tested with initial sawcut lengths 
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of 50.8, 88.9 and 140 mm. Experimental data on applied stress, S, against measured crack 
extension is shown in Figure 15. The critical tearing angle \|/ c , and the sawcut angle \\f{ , 

required in the analysis of the MSD panels were determined from an analysis of these 
specimens and tests on sawcuts reported in reference 10. From the STAGS analyses, a 
sawcut angle of 8.2 degrees and a critical tearing angle of 4.8 degrees were determined to 
best fit these data. 

From the MSD test [16], the only information reported was the applied stress, S, 
against crack-opening displacement at the centerline of the lead crack (crack number 3 in 

Fig. 14). From these results, three peak stresses. Si, S 2 and S 3 , were identified as the first 
crack link up, second crack link up and final failure stress, respectively. The maximum 
stress Sf corresponded to the second peak stress. As these peak stresses could not be 

related to any particular crack length, they are represented by horizontal lines in the Figure 
16. All these test results corresponds to a panel where out-of-plane deformations were 
prevented with guide plates. 



s 

Figure 14. Multiple- site damage (MSD) crack configuration. 
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Figure 15. Comparison of measured and calculated stable crack growth 
in 2024-T3 Alclad aluminum alloy specimens. 

In the STAGS analysis of the MSD panel configuration shown in Figure 14, the 
sawcut and tearing angles determined in the previous section were used. All sawcut tips 
had the same sawcut and tearing angles. The STAGS analysis was able to predict the 
three peak stresses quite accurately, as shown in Figure 16.. For comparison, the results 
from the analysis of a panel with only a lead crack was also shown. The presence of the 
MSD cracks reduced the residual strength of the panel by 30 percent. When the STAGS 
analysis was repeated but allowing out-of-plane deformations (buckling), the stress levels 
at which the cracks link up due to ligament breakage decreased and there was a further 
reduction (10 percent) in the residual strength of the MSD panel. These analyses indicate 
that small MSD cracks combined with buckling can significantly reduce the residual 
strength of MSD panels. 
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Figure 16. Measured and predicted crack link-up for an MSD panel. 

CONCLUDING REMARKS 

Results from the analyses of buckling and stable tearing in thin-sheet materials were 
presented and discussed. Finite-element fracture simulations using the STAGS shell code 
with the critical crack-tip-opening angle (CTO A) fracture criterion produced quite 
accurate results when compared with experimental measurements of crack extension and 
failure loads. The numerical analyses indicated that sheet thickness (B), crack-length-to- 

width (cj/w) ratio, crack configuration, and material fracture toughness had a strong 

influence on the buckling response. The importance of sawcut simulation in predicting 
the residual strength of the panels with single sawcuts and multiple-site-damage cracking 
was highlighted. The STAGS analyses predicted a significant influence of multiple-site 
damage (MSD) cracking on residual strength. There was a further reduction in residual 
strength when the MSD panels were allowed to buckle. Finally, the STAGS code and the 
CTOA fracture criterion may serve as an analysis tool for predicting the residual strength 
of built-up structures subjected to complex loading conditions. 
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Table 1. Material stress-strain properties. 


2024-T3 (LT) 
E = 71.5 GPa 


v = 0.3 
£ 

cr, MPa 

0.00483 

345 

0.015 

390 

0.04 

432 

0.1 

473 

0.16 

491 


7075-T6 


E = 

71.5 GPa 

v = 

0.3 

£ 

a, MPa 

0.00676 484 

0.01 

518 

0.02 

539 

0.04 

560 

0.08 

574 

0.14 

587 

0.2 

588 


2219-T87 
E = 76 GPa 
v = 0.3 


£ 

a, MPa 

0.00455 

346 

0.0074 

378 

0.0183 

411 

0.0423 

453 

0.066 

464 

0.12 

464 


AM355CRT 
E = 207 GPa 
v = 0.3 

£ a, MPa 
0.0074 1534 

0.18 1666 

0.25 1693 


Table 2. Critical crack-tip-opening angles for different materials. 


Material 

Thickness 

Tearing CTOA, \|/ c 

Sawcut CTOA, \|/j 


B, mm 

degrees 

degrees 

2024-T3 (LT) 

2.3 

4.6 

15.0 

2024-T3 Alclad (TL) 

1.0 

4.8 

8.2 

7075-T6 

1.6 

1.9 


2219-T87 

2.5 

3.6 


AM355CRT 

0.5 

4.4 
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